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Based on the quark coalescence model, we derive relations among the momentum anisotropies 
of mesons and baryons in relativistic heavy ion collisions from a given, but arbitrary azimuthal 
distribution for the partons. Besides the familiar even Fourier coefficients such as the elliptic flow, 
we also pay attention to odd Fourier coefficients such as the directed flow, which has been observed 
at finite rapidity even at RHIC energies. 
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I. INTRODUCTION 

The systematic analysis of anisotropies of particle mo- 
mentum distributions in the plane perpendicular to the 
beam direction in non-central collisions of heavy ions of- 
fers the possibility to understand both the dynamics of 
the collisions and the properties of the initially produced 
hot and dense matter [1] . In the low transverse momen- 
tum region, this momentum anisotropy is a result of the 
hydrodynamic expansion that is driven by anisotropic 
pressure gradients in the system [2,3]. Particles of in- 
termediate to high energies do not follow this collective 
dynamics, but instead radiate more or less energy, de- 
pending on how much matter they traverse on their way 
out of the fireball. As a result, these particles also ac- 
quire a strong azimuthal anisotropy in their momentum 
distributions [4]. For heavy-ion collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC), it has been shown in 
hydrodynamical models [5-8] that the lower-order mo- 
mentum anisotropy, such as the elliptic flow, is sensitive 
to the equation of state of the quark-gluon plasma formed 
during the initial stage of the collision as well as the tran- 
sition region to hadronic matter in the later stages. In 
transport models [9-11], the hadron elliptic flow is shown 
to depend on the parton scattering cross sections in the 
initial partonic matter. For hadrons with high transverse 
momentum, the elliptic anisotropy can further provide in- 
formation on the energy density of the initial hot matter 

Although the hadron anisotropies become smaller with 
increasing order, they seem to be sensitive to the initial 
collision dynamics as well [13-15]. Recent experimental 
results from the STAR Collaboration [16] have demon- 
strated that some of the higher-order harmonics are still 
measurable. The experimental data further indicate that 
there is a scaling relation among the hadronic anisotropy 

coefficients v n , i.e., v n (pr) ~ I^u^^Cpt)- To gain in- 
sight in the origin of this scaling relation, we relate in the 
present paper the hadron anisotropies to those of quarks 
using the quark coalescence model. In particular, a naive 
quark coalescence model [17,18], which ignores the mo- 



mentum distribution of quarks inside hadrons, will be 
used. Recently, this model has been shown to give a 
natural explanation for another scaling law observed in 
experiments: the momentum dependence of the elliptic 
anisotropy of identified particles (with the exception of 
pions) becomes similar if both anisotropy and momen- 
tum are divided by the number of constituent quarks in 
the hadron, i.e., two for mesons and three for baryons. 
In Ref . [22] , it was shown that the deviation of the pion 
elliptic flow from the scaling law could be explained by 
taking into account the effect of resonance decays as well 
as the momentum distribution of quarks in hadrons as in 
the more realistic quark coalescence model [19-21]. 



II. ANISOTROPIES IN PARTON AND HADRON 
SPECTRA 

The azimuthal distribution of the particle transverse 
momentum spectra (or the probability of having a parti- 
cle with a given transverse momentum px and azimuthal 
angle tp) can in general be expressed in terms of the 
Fourier series 



f(<p,Pr) oc 1 + 2 22v n (p T ) cosnip 



(2.1) 



where v n denotes the n-th order momentum anisotropy. 
In the above, we have chosen the beam direction for the pr- 
axis and the impact parameter for the x-axis. Sine terms 
which are in general possible in a Fourier expansion al- 
ways vanish in these coordinates because of the reflection 
symmetry with respect to the reaction plane (the (x, z) 
plane). For collisions of equal nuclei, the odd Fourier 
coefficients vanish at mid-rapidity due to the additional 
symmetry with respect to the y-axis. This symmetry 
is lost at forward and backward rapidities as a result of 
spectator matter from the initial nuclei, which eventually 
leads to a tilt of the reaction zone out of the transverse 
plane. Collisions of non-identical nuclei such as S+Au 
break this symmetry even at mid-rapidity. It is there- 
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fore of great interest to keep the odd coefficients in our 
general considerations. 

In the simplest coalescence model [17], mesons of trans- 
verse momentum px are formed by quarks and antiquarks 
of half the transverse momentum pr/2 traveling in the 
same direction. Since the azimuthal distribution func- 
tion Eq.(2.1), after normalizing appropriately, can be in- 
terpreted as a probability distribution in the azimuthal 
angle, the distribution function for mesons is therefore 
given by the square of that for partons, and it can also 
be expressed in terms of a Fourier series, i.e., 



F((p) oc f((p) 2 oc 1 + 2 V n cos n<p , 



(2.2) 



n=l 



where f{ip) given by Eq.(2.1) is taken to be the azimuthal 
distribution of parton transverse momentum spectra, and 
V n denotes the n-th order meson anisotropy. Following 
the calculations outlined in the Appendix, we arrive at 
the following expression for V n at momentum px in terms 
of the parton anisotropy v n at momentum px/2: 



baryon anisotropy V n } s on the left side of the equation 
are evaluated. 

Dominant and clearly visible in the expressions above 
is the linear connection between the n-th harmonic of the 
partonic sector and the n-th harmonic of the hadron sec- 
tor, scaled by a factor 2 and 3 for mesons and baryons, 
respectively. Corrections to these linear relations are at 
least of the order v„ and should therefore be small. This 
gives us a first hint that if the experimental data ex- 
hibit a strong V4 component in the particle momentum 
anisotropy as was predicted in Refs. [15,23] and recently 
confirmed by the STAR collaboration [16], it is then a 
true physical component in the transverse distribution, 
and not just an artifact or trivial byproduct of the ex- 
isting large V 2 signal. We will return to this discussion 
repeatedly. 

We now proceed to study Eqs. (2.3) and (2.6) under 
the assumption that higher-order coefficients are small 
corrections in the distribution and can thus be neglected 
from the sums at certain order in n. 



1 ( n ~ ^ 

V n = — I 2 V n + ^2 V i V n-i + 2 ^2 V * V n+t 
\ t=l i=l 



with 



JV = l + 2^ 2 . 



(2.3) 



(2.4) 



Similarly, we obtain the Fourier spectra of baryons by 
taking the azimuthal distribution of the parton trans- 
verse momentum spectra to the third power, i.e., 



F(ip) oc f(ip) 3 oc 1 + 2 V n cos nip . 



(2.5) 



n=l 



The calculations given in the Appendix then lead to a 
somewhat more complicated expression for the n-th order 
baryon anisotropy V n , 



n-l 



K = i(3w„+3^ V i V n ^ l + 6 V i V n+i 
\ i—l i—l 

oo oc oo n-\-i — 1 

i—l j — 1 i—l j—1 



n—1 n—i—1 
+ X X V ^ V 3 V n-(t+j) 

i=i j=i 



(2.6) 



with 



N = 1 + 6 «? + 6 Yl X • (2-7) 

i—l i—l j=l 

In the above, the parton anisotropy w„'s on the right 
hand side of the equation are determined at a momen- 
tum that is 1/3 of the baryon momentum at which the 



III. ANISOTROPIES AT MID RAPIDITY 

In the first two years of RHIC physics, a vast amount 
of data from Au+Au collisions was collected primarily 
at midrapidity. Due to the symmetry of those systems, 
the odd Fourier coefficients in the expansion of the az- 
imuthal distribution of particle transverse momentum 
vanish both on the partonic and hadronic sector. 

Meson anisotropy at midrapidity: Keeping terms 
up to the order vg and neglecting higher-order correc- 
tions, we find from Eq. (2.3) that for mesons (up to 
order Vq) 



V 2 = ^-( 2 «2 + 2w 2 w 4 + 2w 4 u 6 ) , 
V 4 = ^(2v 4 + vl + 2v 2 v 6 ), 



V 6 



N 



(2v 6 + 2v 2 v A ) . 



(3.1) 
(3.2) 
(3.3) 



with N = 1 + 2{v\ + v\ + vl). 

Interestingly, we notice that a V4 of the mesons can be 
generated out of the partonic sector even if the partonic 
distribution has only an elliptic component v 2 , but no w 4 
contribution. This is an effect of the non-linear way how 
coalescence fuses partons into hadrons (see Eq. (2.2)). 
The resulting V4 is however small, namely for a maximum 
parton anisotropy v 2 of about 10% as required to describe 
the experimental hadron anisotropy one obtains a V4 ~ 
1%. 

Neglecting products of 3 or more w„'s in Eq.(3.1), we 
find 



1 v 4 
24 



(3.4) 
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Experimentally [16], this ratio (measured for all charged 
hadrons) is approximately constant as a function of trans- 
verse momentum, and has a value of about 1.2. There- 
fore, we can be quite certain that there are higher-order 
harmonics in the Fourier expansion of the azimuthal de- 
pendence of the parton transverse momentum distribu- 
tion. In fact we can deduce from this observation that 
va{pt) ~ 2w|(pt)- A finite parton V4 is indeed seen in 
the predicted parton anisotropics [23] from the AMPT 
(A Multi-Phase Transport) model [24], although it scales 
with v 2 with a scaling factor of only about 1 . 

Baryon anisotropy at midrapidity: For the 
Fourier coefficients of the baryon anisotropy, we again 
consider coefficients up to 6th order. From Eq. (2.6), we 
obtain 



V 2 = ^(3^2 + 6v 2 v 4 + 6u 4 « 6 + 



(3.5) 



+3w| + 3w|w 6 + 6v 2 vl + 6v 2 v% + 3wf i> 6 ) , 



-6v 2 V4 + 6v 2 V4Vq + 3v 4 + 6v 4 v 6 ) 



V 6 = ^"(3^6 + 6v 2 v 4 



(3.6) 



(3.7) 



withiV= l + 6(v%+vl + vl) + 6(v%V4 + 2v2V 4 v 6 ). We note 
that to obtain the anisotropy values at momentum pt on 
the left hand side for baryons , we have to use the parton 
anisotropy coefficients at pr/3 on the right hand side. As 
already discussed before, there is a linear connection of 
quark and baryon anisotropics [17], Vn(pr) ~ 3v n (pT/3) 
with corrections which are of the order v„ . 

Again, we express the ratio of the 4-th order coefficient 
V4 to the square of the elliptic coefficient V 2 up to second 
order in partonic v n , and the result is 



1 v 4 



(3. 



This confirms the earlier argument, that the value ob- 
served by the STAR collaboration can only be achieved 
if there is a significant V4 component on the quark sector. 

Together, Eqs. (3.4) and (3.8) have the immediate 
consequence that independent of the detailed ccntral- 
ity or momentum dependence of the quark anisotropics, 
there is a general relation between baryon and meson 
anisotropies, 



(3.9) 



An experimental investigation of the ratios for identified 
particles thus proves promising to further approve or dis- 
prove the quark coalescence picture, which we leave here 
as a challenge for future experimental analysis and cor- 
respondingly for the coalescence model. 



IV. ANISOTROPIES FOR NON-ZERO RAPIDITY 
OR NON-IDENTICAL COLLISION PARTNERS 

At non-zero rapidity the symmetry x <-> —x, which is 
apparent for collisions of equal nuclei, is broken. This 
allows for the presence of odd Fourier coefficients both in 
the partonic and hadronic sectors. This specific symme- 
try is also lost at midrapidity for collisions of non-equal 
nuclei, such as S+Au. Even if one does not expect to cre- 
ate a thermalized system in smaller systems, such as in 
last years d+Au run at RHIC, the following discussion is 
also useful to quantify and relate the different azimuthal 
coefficients at midrapidity. 

Meson anisotropy: We restrict our discussion to the 
first four coefficients v\ to v 4 in the partonic distribution. 
It then follows directly from Eq.(2.3) that 



Vi = (2vi + 2viv 2 + 2v 2 v 3 + 2v 3 v 4 ) 

l\y 

V 2 = -t- (2v 2 + v\ + 2v x v z + 2v 2 V4) , 

IV71 



V3 = — (2v 3 + 2v x v 2 + 2wi« 4 ) , 

l\y 



v 4 



N, 



(2v 4 +vl + 2v x v z ) , 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



vl). 



with N y = 1 

As for the even coefficients in the section before, we see 
that the mesonic V n is linearly related to the partonic v n 
with corrections that are at least quadratic in v n . Even 
more interesting is that here a meson V 2 can be generated 
purely through the partonic directed flow V\ without any 
dipole deformation in the initial matter. 

Another interesting point is the relation between har- 
monics of different orders. As we have already discussed, 
the experimental data at midrapidity indicate a scaling of 
V4 with V 2 as also insinuated by the coalescence model. 
At forward rapidity we can investigate a similar relation 
for V3. We notice that a relation between V3 and is 
unlikely because the product of odd harmonics gives rise 
to an even one, and this is indeed the reason why there 
are no v 2n+1 in V 2n +i while they are present in V 2n . This 
suggests a possible relation between V3 and the odd-even 
combination V\V 2 . From the above formulas, we find that 



V1V2 



1 v 3 



2 v\v 2 



1 V4 
~2V 2 



(4.5) 



The presence of non-negligible odd harmonics at for- 
ward rapidity can also reflect on the ratio of fourth to 
second component by adding a new term to Eq.(3.4), 



i.e., 



V4 1 1 V 4 1 V\Vz 

Vi ~ 4 + 2^| + 2^|" 



(4.6) 



Baryon anisotropy: The expressions for the baryon 
anisotropy coefficients are rather lengthy but become eas- 
ily accessible when dropping the terms that contain prod- 
ucts of three v n coefficients, i.e., by dropping the last two 
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lines of Eq.(2.6). The resulting approximate expressions 
are then 



AS, 



and 



N,. 



[3«i 


+ 6(viv 2 + v 2 v 3 + V3V4)], 


(4.7) 


[3u 2 


+ 3{vl + 2v 1 v 3 + 2v 2 v 4 )}, 


(4.8) 


[3u 3 


+ 3{2v 1 v 2 + 2v 1 v 4 )}, 


(4.9) 


[3w 4 


+ 3(2«i«3+^)], 


(4.10) 


1 + 


&{v{ + vl + vl + vl). 


(4.11) 



Again, we notice that the n-th coefficient of the baryon 
spectra is linearly related to the n-th coefficient in the 
partonic sector, and higher-order terms are corrections 
of the order However, the above expressions also 
show that even harmonics in the baryon spectra can be 
generated from only odd coefficients in the quark sector, 
even in the absence of even partonic harmonics. 

As in the case of mesons, it is of interest to look at the 
ratio between the third and the first two components in 
the azimuthal distribution of baryon momentum spectra. 
In the coalescence picture, this is given by 



V 3 



1 V3 f + ^a± 

3 viv 2 \ 3v 2 



(4.12) 



This relation allows us also to investigate the presence of 
a non- vanishing v 3 at quark level. In fact, if v 3 is equal 
to zero the ratio V3/V1 V 2 — 0.7. Therefore, the existence 
of a quark 113 component can be investigated by looking 
at the magnitude of this ratio from experimental data, 
as done for V4. 

Similar to the case of hadrons at midrapidity in col- 
lisions of equal mass nuclei, where baryon V^/V 2 and 
meson Vi/V 2 are related by Eq.(3.9), we have for the 
case of non- vanishing odd-order anisotropies the follow- 
ing relation: 

if the small U4 contribution is neglected in Eqs.(4.5) and 
(4.12). 

We want to point out that the above considerations 
for both mesons and baryons are applicable in a rapidity 
region where the hadrons are produced by coalescence 
from the previously established deconfined phase. These 
formulas therefore do not capture the region of highest ra- 
pidity, the fragmentation region, which is located around 
beam rapidity. 



V. SUMMARY 

Using a naive quark coalescence model for hadron pro- 
duction from the quark-gluon plasma formed in rela- 
tivistic heavy ion collisions, which only allows partons 
of equal momenta to recombine, we have expressed the 
Fourier coefficients of the azimuthal distribution of me- 
son and baryon momentum spectra in terms of those of 
partons. Neglecting higher-order terms, simple relations 
are found between the ratio of the hadrons 4th order 
term to the square of their 2nd order to the correspond- 
ing ratio for the partons. Furthermore, a simple relation 
is obtained between such ratios for baryons and mesons, 
i.e., Eq. (3.9). This poses an important prediction of the 
coalescence model independently of other effects such as 
flow or jet-quenching. For hadrons at finite rapidities or 
spectra from collisions of non-identical nuclei, odd-order 
terms are nonzero. We have shown that the ratio of their 
third-order term to the product of their lst-order and 
2nd-order terms is also simply related to the correspond- 
ing ratio for the partons. To account for the observed 
large ratio (~ 1.2) of hadron 4th-order to the square of 
2nd-ordcr anisotropy, a larger value (~ 2) for the corre- 
sponding ratio for the partonic anisotropy is needed. Al- 
though this value is larger than that seen in the AMPT 
model (~ 1) [23], which is based on a more realistic coa- 
lescence picture that includes both the quark momentum 
distributions inside hadrons and also the effect of reso- 
nance decays, the present schematic study allows one to 
have a simple understanding of the effect of higher-order 
partonic anisotropies on those of the observed hadrons. 



APPENDIX 

To derive the harmonic coefficients on the hadronic sec- 
tor, consider an arbitrary distribution function f(<p) on 
the interval ip £ [0, 2ir] expressed in terms of its Fourier 
series 



/(v) = fk cos k(p . 



(5.1) 



k=0 



Squaring this function, as it is required for the meson 
distributions, we obtain 



F (<fi) = f(<f) 2 = ^2^2.fkficosk(pcosl(p 

k=0 1=0 

1 ' * 

= o f k h [ cos ( fc +0^ + cos(fc-Z)v] • 



k=0 1=0 



From this we then find the term independent of the az- 
imuthal angle 



(5.2) 



k=0 
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and the amplitude of the nth harmonic 

^ n oo 
fe=0 fe=0 

Substituting / = 1 and = 2i^, we obtain the expres- 
sions given in Eq.(2.3). 

For baryons, we have to apply the third power to the 
function f(<p), which leads to 

F{<p) = f{pf (5.4) 

oo oo oc 

= XX X fkflfm cos kp cos Ip cos mip 

k=0 1=0 m=0 

= \ X fkflfm [cos(k+l + m)p + cos(fc+Z- m)<p 

k,l,m 

+ cos(k — l + m)p + cos(k — l — m)p} 
= \ X fkfifm[cos(k+l+m)<p> + 3cos(k+l-m)ip] , 

where the last equality follows from the repeatedly oc- 
curring summation indices in the line above. Again, we 
can quickly separate the terms with a specific amplitude 
and find 

1 3 °° 

k,l 

and 

^ oo oo 3 °° nJrk 

Fn = J X X fkflfn+k+l + 7 X X fkflfn+k-l 
k=0 1=0 k=0 1=0 

^ n n—k 

+ iEE fkflfn-k-l ■ (5.6) 
fe=0 (=0 

With the substitution / = 1 and fk — 2vk, one then 
obtains Eq.(2.6). 
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